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In this paper we investigate coherent and squeezed quantum states of phonons. The latter allow 
the possibility of modulating the quantum fluctuations of atomic displacements below the zero- 
point quantum noise level of coherent states. The expectation values and quantum fluctuations of 
both the atomic displacement and the lattice amplitude operators are calculated in these states — in 
some cases analytically. We also study the possibility of squeezing quantum noise in the atomic 
displacement using a polariton-based approach. 
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I. INTRODUCTION 

Classical phonon opticsEI has succeeded in producing 
many acoustic analogs of classical optics, such as phonon 
mirrors, phonon lenses, phonon filters, and even "phonon 
microscopes" that can generate acoustic pictures with a 
resolution comparable to that of visible light microscopy. 
Most phonon optics experiments use heat pulses or super- 
conducting transducers to generate incoherent phonons, 
which propagate ballistically in the crystal. These ballis- 
tic incoherent phonons can then be manipulated by the 
above-mentioned devices, just like in geometric optics. 

Phonons can also be excited phase- coherently. For in- 
stance, coherent acoustic waves with frequencies of up 
to 10 10 Hz can be generated by piezoelectric oscillators. 
Lasers have also been used to generate coherent acoustic 
and optical phonons through stimulated Brillouin and 
Raman scattering experiments. Furthermore, in recent 
years, it has bees, possible to track the phases of coherent 
optical phononsB, due to the availability of femtosecond- 
pulse ultrafast lasers (with a pulse duration shorter than 
a phonon period)0, and techniques that can measure op- 
tical reflectivity with accuracy of one part in 10 6 . 

In most situations involving phonons, a classical de- 
scription is adequate. However, at low enough temper- 
atures, quantum fluctuations become dominant. For ex- 
ample, a recent studyu shows that quantum fluctuations 
in the atomic positions can indeed influence observable 
quantities (e.g., the Raman line shape) even when tem- 
peratures are not very low. With these facts in mind, and 
prompted by the many exciting developments in classi- 
cal phonon optics, coherent phonon experiments, and (on 
the other hand) squeezed states of lightO, we would like to 
explore phonon analogs of quantum optics. In particular, 
we study the dynamical and quantum fluctuation prop- 
erties of the atomic displacements, in analogy with the 
modulation of quantum noise in light. Specifically, we 
study single-mode and two-mode phonon coherent and 
squeezed states, and then focus on a polariton-based ap- 
proach to achieve smaller quantum noise than the zero- 



point fluctuations of the atomic lattice. 

The concepts of coherent and squeezed states were 
both originally proposed in the context of quantum op- 
tics. A coherent state is a phase-coherent sum of number 
states. In it, the quantum fluctuations in any pair of 
conjugate variables are at the lower limit of the Heisen- 
berg uncertainty principle. In other words, a coher- 
ent state is as "quiet" as the vacuum state. Squeezed 
statesB are interesting because they can have smaller 
quantum noise than the vacuum state in one of the con- 
jugate variables, thus having a promising future in dif- 
ferent applications ranging from gravitational wave de- 
tection to optical communications. In addition, squeezed 
states form an exciting group of states and can provide 
unique insight into quantum mechanical fluctuations. In- 
deed, squeezed states are now being explored in a va- 
riety of non-quantum-optics systems, including classical 
squeezed statesO. 

In Sec. II we introduce some quantities of interest 
and study the fluctuation properties of the phonon vac- 
uum and number states. In Sees. Ill and IV we inves- 
tigate phonon coherent and squeezed states. In Sec. V 
we propose a novel way of squeezing quantum noise in 
the atomic displacement operator using a polariton-based 
mechanism. The Appendix summarizes the derivation of 
the time evolution of the relevant operators in this polari- 
ton approach. Finally, Sec. VI presents some concluding 
remarks. 



II. PHONON OPERATORS AND THE PHONON 
VACUUM AND NUMBER STATES 

A phonon with quasimomentum p = hq and branch 
subscript A has energy e q ,\ = hujq\; the correspond- 
ing creation and annihilation operators satisfy the boson 
commutation relations 

[bq'X', 6q A ] = <W^AA' ' [ fe qA> Vv] = . (I) 

The atomic displacements Ui a of a crystal lattice are 
given by 
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Here refers to the equilibrium lattice positions, a to 
a particular direction, and Q q is the normal-mode am- 
plitude operator 



Q x /JL(, 
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An experimentally observable quantity is the real part of 
the Fourier transform of the atomic displacement: 



Re[w Q (q)] = E 



8m w qA 

+^qa(6-qA + 0} 



{^JfcqA + ^qJ 



(4) 



For simplicity, hereafter we will drop the branch sub- 
script A, assume that £/ qa is real, and define a q-mode 
dimcnsionlcss lattice amplitude operator: 



u(±q) 



6 q + bi^ + 6_ q 



(5) 



This operator contains essential information on the lat- 
tice dynamics, including quantum fluctuations. It is the 
phonon analog of the electric field in the photon case. 

Let us first consider the phonon vacuum state. When 
no phonon is excited, the crystal lattice is in the phonon 
vacuum state |0). The expectation values of the atomic 
displacement and the lattice amplitude are zero, but the 
fluctuations will be finite: 
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Let us now consider the phonon number states. The 
eigenstates of the harmonic phonon Hamiltonian are 
number states which satisfy 6 q |". q ) = y^ql^q — !)■ The 
phonon number and the phase of atomic vibrations are 
conjugate variables. Thus, due to the uncertainty prin- 
ciple, the phase is arbitrary when the phonon number 
is certain, as it is the case with any number state |n q ). 
Thus, in a number state, the expectation values of the 
atomic displacement (n q |ui Q |n q ) and q-mode lattice am- 
plitude (n q |u(±q)|n q ) vanish due to the randomness in 
the phase of the atomic displacements. The fluctuations 
in a number state \n a ) are 



III. PHONON COHERENT STATES 

A single-mode (q) phonon coherent state is an eigen- 
statc of a phonon annihilation operator: 



/3ql/3q) 



(11) 



It can also be generated by applying a phonon displace- 
ment operator _D q (/3 q ) to the phonon vacuum state 



l/3q) 



Dq(/3 q )|0) =exp(/3 q 6 q -/5 q 6 q )|0) 

E Mq 



exp 
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Thus it can be seen that a phonon coherent state is a 
phase-coherent superposition of number states. More- 
over, coherent states are a set of minimum- uncertainty 
states which are as noiseless as the vacuum state. Co- 
herent states are also the set of quantum slates that best 
describe the classical harmonic oscillatorso. 

A single-mode phonon coherent state can be generated 
by the Hamiltonian 

H = ftw q (b% + ~) + A q (t) 6 q + A q (i) 6 q (14) 

and an appropriate initial state. Here A q (t) represents 
the interaction strength between phonons and the exter- 
nal source. More specifically, if the initial state is a vac- 
uum state, \ip(0)) = |0), then the state vector becomes a 
single-mode coherent state thereafter 



where 



|^)) = |A q (i)e- 



A q (i) = -~/ A q (r)e— dr 



(15) 
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is the coherent amplitude of mode q. If the initial state 
is a single-mode coherent state IV'(O)) = |a q ), then the 
state vector at time t takes the form 



|V(*)) = |{A q (t)+o q }e- fa "*) 



(17) 



which is still a coherent state. 

In a single-mode (q) coherent state |A q (i) e _l " q *), 
(ui a (t)}coh and (u(±q)) co h are sinusoidal functions of 
time. The fluctuation in the atomic displacements is 



JV 



2Nm cu c 



(18) 



<(Au JQ ) 2 ) r 
((A^(±q)) 2 



h\U, 



N 

qa | X ^ 



h\Ua> a \ 



Nra Ld n 



2NmuJn' c 

q'#q 



1 num — 2 -(- 27^ q . 



(9) 
(10) 



The unexcited modes are in the vacuum state and thus 
all contribute to the noise in the form of zero point fluc- 
tuations. Furthermore, 



((A U (±q)) 2 ) coh = 2. 



(19) 
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From the expressions of the noise ((Awio,) 2 ) co h and 
((Au[±q)] 2 ) co h, it is impossible to know which state (if 
any) has been excited, while this information is clearly 
present in the expression of the expectation value of 
the lattice amplitude (it(±q)) co h- These results can 
be straightforwardly generalized to multi-mode coherent 
states. 



IV. PHONON SQUEEZED STATES 

In order to reduce quantum noise to a level below the 
zero-point fluctuation level, we need to consider phonon 
squeezed states. Quadrature squeezed states are gener- 
alized coherent states!! Here "quadrature" refers to the 
dimensionless coordinate and momentum. Compared to 
coherent states, squeezed ones can achieve smaller vari- 
ances for one of the quadratures during certain time in- 
tervals and are therefore helpful for decreasing quantum 
noise. Figures [l] and || schematically illustrate several 
types of phonon states, including vacuum, number, co- 
herent, and squeezed states. These figures are the phonon 
analogs jpf the illuminating schematic diagrams used for 
photons! . 

A single-mode quadrature phonon squeezed state is 
generated from a vacuum state as 

|a q ,0 = Aj(oO£ q (0|0>; (2°) 

a two-mode quadrature phonon squeezed state is gener- 
ated as 

l^qi i Q; q2 ' 

= A u (a qi )D q2 (a q2 )S qi , q2 (O|0) ■ (21) 

Here Z? q (a q ) is the coherent state displacement operator 
with a q = |a q |e I( ^, 

S q (£)=exp(^6 q -i& q 2 ) , (22) 

S qi ,q 2 (0 = cxp (r& q A 2 - ^L&L) , (23) 

are the single- and two-mode squeezing operator, and 
£ = re 10 is the complex squeezing factor with r > and 
< 9 < 2tt. The squeezing operator 5 qiiq2 (^) can be 
produced by the following Hamiltonian: 

#qi,qa = ^qi^q^qi + ^q 2 & q 2 & q2 

+c(*KX 2 + n < )Mq 2 - (24) 

The time-evolution operator has the form 

U(t) = exp (~S t\ exp [f (t)6 qi 6 q2 - fflb^b^] , 

(25) 

where 

Ho = hcjq, 6 qi fe qi + hu;^ b\ 2 6 q2 , (26) 



and 

£(i) = i/ C(r)e l( ^ + ^ )T dT. (27) 

ft J-oo 

Here is the squeezing factor and ((t) is the strength 
of the interaction between the phonon system and the ex- 
ternal source; this interaction allows the generation and 
absorption of two phonons at a time. The two-mode 
phonon quadrature operators have the form 

A(q, -q) = 2- 3 / 2 (& q + 6 q + 6_ q + 6^) (28) 
= 2~ 3 / 2 M (±q) , (29) 
P(q, -q) = -i2~ 3 / 2 (& q - b{ + &_ q - 6l q ) . (30) 

We have considered two cases where squeezed states 
were involved in modes ±q. In the first case, the system 
is in a two-mode (±q) squeezed state |a q ,a_ q ,£), (£ = 
re t6 ), and its fluctuation is 

([Au(±q)] 2 ) = 2 (V 2 - cos 2 °- + e 2r sin 2 . (31) 

In the second case, the system is in a single-mode 
squeezed state |a q ,£) (a q = |a q |e l( ^) in the first mode 
and an arbitrary coherent state |/3_ q ) in the second mode. 
The fluctuation is now 

([A U (±q)] 2 ) = l + e 2r sin 2 U+j) 

+e- 2r cos 2 U+^j ■ (32) 

In both of these cases, ([Au(±q)] 2 ) can be smaller than 
in coherent states (see Fig. 2). 

V. POLARITON APPROACH 

Phonon squeezed states can be generated through 
phonon-phonon interactions. This will be discussed 
elsewherea. Here we focus on how to squeeze quan- 
tum noise in the atomic displacements through phonon- 
photon interactions. When an ionic crystal is illuminated 
by light, there can be a strong coupling between photons 
and the local polarization of the crystal in the form of 
phonons. Photons and phoncps with the same wave vec- 
tor can thus form polaritonslllJ. Although now phonons 
and photons are not separable in a polariton, we can still 
study the quantum noise in the atomic—displacements. 
Let us consider the simplest Hamiltoniar£3 describing the 
above scenario: 

Spolaxiton = ^ {-Elk a\fl\s. + -E^k&k &k 
k 

+E 3k (albk - ak& k - a k fo-k + a- k & k) } > (33) 
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where 



VI. CONCLUSIONS 



Eik = hck , 



E 3 k = i 



. f h 2 ckuJoX 



1/2 



(34) 
(35) 

(36) 



Here k is the wave vector for both photons and phonons 
and Wo is the bare phonon frequency. \ is the dimension- 
less dielectric susceptibility of the crystal (the strength 
of the phonon-photon interaction) defined by 



XWo£oE = p + ugp . 



(37) 



where E is the electric field of the incoming light and 
P is the polarization generated by optical phonons in 
the crystal. In -f/poiariton, the two free oscillator sums 
correspond to free photons and free phonons, while the 
mixing terms come from the interaction E • P between 
photons and phonons. The phonon energy i?2k has been 
corrected as u>o is substituted by u>q^/1 + \, so that we 
have "dressed" phonons. 

Our goal is to compute the fluctuations of the lattice 
amplitude operator u(±k, t) = bk(t) + b_ k (t) + 6_k(i) + 
ln a two-mode (±k) coherent state |a q , a_ q ), its 



variance is ([Au(±q)] z 



' coh 



2. Therefore, if at any given 



time we obtain a value less than 2, the lattice amplitude 
of the relevant mode is squeezed. In our calculation, we 
diagonalize the polariton Hamiltonian and find the time- 
dependence of u(±q). The Appendix presents in more 
detail the derivation of the time-evolution of u(±q). 

Our results show that the fluctuation property of 
w(±q) sensitively depends on the t = initial state \ip{0)} 
of both phonons and photons. Our results are summa- 
rized in Table || and some numerical examples are shown 
in Fig. |j| These calculations focus on the case where ck 
is close to wo (the bare phonon frequency, which is typi- 
cally ~ 10 THz for optical phonons) and thus our typical 
time is ~ 0.1 ps. More specifically, squeezing effects in 
w(±k) are relatively strong for either one of the following 
two sets of t = initial states: (i) photon and phonon co- 
herent states, or (ii) single-mode photon squeezed state 
and phonon vacuum state. For instance, the maximum 
squeezing exponent r is 0.015 when the incident pho- 
ton state has a squeezing factor £ = 0.1e 2lcfc * (where ck 
is the photon frequency). On the other hand, with an 
initial two-mode photon (±k) squeezed state and two- 
mode (±k) phonon vacuum state, the squeezing effect in 
w(±k) is weak. We have also used initial conditions with 
a single-mode photon squeezed state and thermal states 
in the twophonon modes. 

Figure |j shows the temperature dependence of the 
squeezing effect for several values of the dielectric suscep- 
tibility x of the crystal. Our numerical results show that 
squeezing effects are quickly overshadowed by the ther- 
mal noise for small x, while for large x ( e -g-? X = 0-5) the 
squeezing effect can exist up to T w 250 K, as illustrated 
in Fig. | 



In conclusion, we have investigated the dynamics and 
quantum fluctuation properties of phonon coherent and 
squeezed states. In particular, we calculate the exper- 
imentally observable time evolution and fluctuation of 
the lattice amplitude operator u(±q). We show that the 
(w(±q)) are sinusoidal functions of time in both coherent 
and squeezed states, but the fluctuation ((Au(±q)) 2 ) in 
a squeezed phonon state is periodically smaller than its 
vacuum or coherent state value 2. Therefore, phonon 
squeezed states are periodically quieter than the vac- 
uum state. In the polariton approach to squeezing, we 
calculate the atomic displacement part of a polariton, 
and prove that the fluctuations of the associated lattice 
amplitude operator can be squeezed for different com- 
binations of initial photon and phonon states and large 
enough (x > 0.1) interaction strength. 

It is difficult to generate squeezed states because they 
have noise levels which are even lower than the one for 
the vacuum state. Indeed, the experimental and theoret- 
ical development of photon coherent and squeezed states 
took decades. Likewise, the experimental realization of 
phonon squeezed states might require years of further 
theoretical and experimental work. Nevertheless, we be- 
lieve that theoretical results in quantum phonon optics 
can help the development of the corresponding experi- 
ments. We hope that our effortH into this very rich prob- 
lem will lead to more theoretical and experimental devel- 
opments in the still unexplored area of quantum phonon 
optics and the manipulation of phonon quantum fluctu- 
ations. 
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APPENDIX A: TIME-EVOLUTION OF THE 
PHONON OPERATORS IN A POLARITON 

To derive the time-evolution of the phonon operators, 
we need to first diagonalize the polariton Hamiltonian 
-ffpoiariton- For this purpose, we introduce the polariton 
operators in terms of the phonon and photon opera- 
tors 6k and ctk 

a ik = w,a k + Xi 6 k + Ui a f _ k + z t &l k , i = l,2. (Al) 
If we write 



a = (aik, a 2 k, a>\ _ k , al _ k ) T . 
a = (a k , 6k, aL k , & _k) T > 



(A2) 
(A3) 
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the above relation can be written in a matrix form 

a = A ■ a , (A4) 
with its inverse a = A -1 • a. Here A is a matrix given 

by 



I W\ X! J/1 Z\ \ 

W 2 X 2 J/2 z 2 

It ak ak __ifc 

2/! X! 

V 2/2 4 W 2 ^2 / 



(A5) 



In the polariton representation, the Hamiltonian has 
the diagonal form 



^polariton = Yl E ^ ( "ik " lk + ^ 

L- L \ 



+3 



(2) 



C*2k "2k + ^ 



(A6) 



(A7) 



The subindiccs i = 1,2 specify the two polariton 
branches, with different dispersion relations E^ and 

(2) 

E^'. The transformation matrix elements wi, Xi, j/i, 2» 
are determined by requiring that the 's satisfy boson 
commutation relations 



«ik , aj k ,] = % <5kk' , [a»k , «jk'] = , 



so that 



[a lk , iJ] = E^ ] a jk , 



(A8) 



(A9) 



which is true if the two different polariton branches are 
independent of each other. 

In the polariton representation, the Hamiltonian 
-^polariton describes two independent harmonic oscilla- 
tors. From the Heisenberg equation 



iE ( k r> t/h 



we obtain 

aik(t) = aik(O) e 
a 2k (t) = a 2k (0) e-^ ; */'< 
or in a more compact form 

a(f) = U a (i)a(0) 



(A10) 

(Alia) 
(Allb) 

(A12) 



Recall that the matrix form of the canonical transforma- 
tion from the photon and phonon operators (a k and 6 k ) 
to the polariton operators (a k ) is a = A ■ a. Thus, at 
time t the photon and phonon operators can be expressed 
as 



a(t) = A" 1 a(t) = A' 1 U«(t)Aa(0) , 



(A13) 



which provides the time evolution of the photon and 
phonon operators. 
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t = photons 


t — phonons 


Squeezed u(±k)? 


CS(±k) 


CS(±k) 


yes (no) if x > (<) 0.1 


SMST(k), VS(-k) 


VS(±k) 


yes (no) if x > (<) 0.1 


SMST(k), VS(-k) 


TS(±k) 


yes if T<T a ( X ) 


TMST(±k) 


VS(±k) 


weak (no) if x > (<) 0.1 



TABLE I. Different combinations of t = initial states 
(modes ±k) for the polariton approach to lattice amplitude 
squeezing and the corresponding effects in the fluctuations of 
the lattice amplitude operator u(±k). Here CS(k), VS(k), 
TS(k), SMST(k), TMST(±k) refer, respectively, to coherent, 
vacuum, thermal, single-mode and two-mode squeezed states 
in the mode inside the parentheses, k or ±k. T s (x) is the 
temperature below which squeezing is obtained (see Fig. 4). 
By squeezing we mean that the quantum noise of the relevant 
variable is below its corresponding vacuum state value. 
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FIG. 1. Schematic diagram of the uncertainty areas 
(shaded) in the generalized coordinate and momentum 
(X(q, — q), P(q, — q)) phase space of (a) the phonon vac- 
uum state, (b) a phonon number state, (c) a phonon coherent 
state, and (d) a phonon squeezed state. Here X(q, — q) and 
P(q, — q) are the two-mode (±q) coordinate and momentum 
operators, which are the direct generalizations of their cor- 
responding single-mode operators. Notice that the phonon 
coherent state has the same uncertainty area as the vacuum 
state, and that both areas are circular, while the squeezed 
state has an elliptical uncertainty area. Therefore, in the 
direction parallel to the 6/2 line, the squeezed state has a 
smaller noise than both the vacuum and coherent states. 



FIG. 4. Temperature dependence of the minimum fluctua- 
tion min{((Aii[±k)] 2 )} in u(±k) using a polariton mechanism 
for squeezing. The phonon frequency is 10 THz. The initial 
states are: single-mode squeezed state in photon mode k, 
vacuum state in photon mode — k, and thermal state in both 
(±k) phonon modes. The squeezing factor is £ = O.le 2 *'. 
Squeezing can exist up to a temperature T s (x)- For example, 
when x = 0.2, squeezing effects vanish when T > 25 K. On 
the other hand, for stronger photon- phonon interaction (e.g., 
X = 0.5), the squeezing effects can be obtained up to T w 250 
K. 



FIG. 2. Schematic diagram of the time evolution of the 
expectation value and the fluctuation of the lattice ampli- 
tude operator u(±q) in different states. Dashed lines repre- 
sent (u(±q)), while the solid lines represent the envelopes 
(u(±q)) ± ^/{[Au(±ci)] 2 ). (a) The phonon vacuum state 
|0), where (u(±q)) = and <[Au(±q)] 2 ) = 2. (b) A 
phonon number state |n q ,n_ q ), where (w(±q)) = and 
([Att(±q)] 2 ) = 2(n q + 7i_ q ) + 2. (c) A single-mode phonon 
coherent state |a q ), where (u(±q)) = 2|a q |cosaj q t (i.e., a q 
is real), and {(Au[±q)] 2 ) = 2. (d) A single-mode phonon 
squeezed state ja q e _l " qt , £(£)}, with the squeezing factor 
£(t) = re" 2<Wt1 ' and r = 1. Here, (u(±q)) = 2|a q |cosw q t, 
and <[Att(±q)] 2 ) = 2(e~ 2r cos 2 w^t + e 2r sin 2 w q t). (e) A sin- 
gle-mode phonon squeezed state, as in (d) ; now the expecta- 
tion value of u is (u(±q)) = 2|a q |sinw q t, (i.e. a q is purely 
imaginary), and the fluctuation ([Au(±q)] 2 ) has the same 
time-dependence as in (d). Notice that the squeezing effect 
now appears at the times when the lattice amplitude (u(±q)) 
reaches its maxima, while in (d) the squeezing effect is present 
at the times when (it(±q)) is close to zero. 



FIG. 3. Calculated ([Au(±k)] 2 ) versus time for different 
combinations of photon and phonon initial states using a po- 
lariton mechanism for lattice amplitude squeezing. Dashed 
(solid) lines correspond to a susceptibility \ = 0.1 (0.4). 
Time is measured in units of 1/cfc, where ck is the free pho- 
ton frequency. These calculations focus on the case where 
ck is close to uio (the bare phonon frequency, which is typ- 
ically ~ 10 THz for optical phonons) and thus our typical 
time is ~ 0.1 ps. The horizontal lines at ([Au(±k)] 2 ) = 2 
correspond to the noise level of coherent states. Thus, any 
time the fluctuation satisfies ([Ait(±k)] 2 ) < 2 (highlighted), 
the state is squeezed. Different combinations of initial states 
were considered, (a) Photon and phonon coherent states, (b) 
Single-mode squeezed state in photon mode k with squeez- 
ing factor £ = 0.1 and a vacuum state in the photon mode 
— k; both phonon modes are in the vacuum state, (c) Same 
combination of states as in (b), but here £ = 0.1e 2lt . 
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